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MINOS COLLABORATION 

We present an analysis designed to search for Lorentz and CPT violations as pre- 
dicted by the SME framework using the charged current neutrino events in the 
MINOS near detector. In particular we develop methods to identify periodic vari- 
ations in the normalized number of charged current neutrino events as a function 
of sidereal phase. To test these methods, we simulated a set of 1,000 experiments 
without Lorentz and CPT violation signals using the standard MINOS Monte 
Carlo. We performed an FFT on each of the simulated experiments to find the 
distribution of powers in the sidereal phase diagram without a signal. We then 
injected a signal of increasing strength into the sidereal neutrino oscillation prob- 
ability until we found a 5cr deviation from the mean in the FFT power spectrum. 
By this method, we can establish upper limits for the Lorentz and CPT violating 
terms in the SME. 



1. Introduction 

At experimentally accessible energies, signals for Lorentz and CPT viola- 
tion can be described by the Standard Model Extension (SME) 1 ' 2 , a theory 
whose foundation is the Standard Model (SM) of particle physics. Since the 
SM is expected to be the low-energy limit of a more fundamental theory 
that unifies quantum physics and gravity at the Planck scale, m p ~ 10 19 
GeV, the violations of Lorentz and CPT symmetries predicted by the SME 
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provide a link to Planck scale physics. Assuming a quantum-gravity origin 
for the violations, however, suggests their magnitude in the accessible en- 
ergy limit is suppressed by a factor of 1CP 17 or more 3 . Despite this huge 
suppression, these low-energy signatures from new physics at the Planck 
scale can be probed with current experimental technologies. 

The SME framework predicts several unconventional neutrino signals, 
among which is one that arises from the dependence of the neutrino os- 
cillation probability on the direction of neutrino propagation 3 . For exper- 
iments like MINOS 4 with both beam neutrino source and detector fixed 
on the Earth's surface, the Earth's sidereal rotation causes the direction 
of neutrino propagation p to change with respect to the Sun-centered in- 
ertial frame in which the SME is formulated 5 . The theory predicts that 
this rotation introduces a sidereal variation in the number of neutrinos de- 
tected from the beam. In this paper we use simulated beam neutrinos in 
the MINOS near detector 6 to develop methods to search for this sidereal 
signal. 

We developed our search algorithms by injecting an SME signal into 
our simulated data set of MINOS near detector neutrinos and then 
developing methods to find it. According to SME, the probability that a 

oscillates from — > v T over a distance L due to Lorentz and CPT 
violation is given by 3 

P v ^ Vt ~ L 2 \(C)^ Vt + (A c )^ Vt cos (Lu m T e ) (1) 
+ {A S ) V ^ T sin (w e T e ) + (-Bc)^ cos (2w®T ffi ) 
+(B a )„„„ T sin(2^ ffi T ffi )] 2 , 

where = 2ir/(2'i h 56 m 04.0982 s ) is the Earth's sidereal frequency and 
T e is the Local Sidereal time of the event. In this equation, the expressions 
for (A C )L, (A S )L, and (C)L include both CPT and Lorentz violating terms; 
the expressions for (B C )L and (B S )L include only Lorentz violating terms. 
There are no terms that depend on CPT violation alone because CPT 
violation implies Lorentz violation 7 . In cq.(l), the CPT violating terms 
depend only on L and the Lorentz violating terms depend on L x E Vfi , 
where E v is the energy of the neutrino. This unconventional behavior is 
to be compared with the L/E Vfi dependence of the oscillation probability 
when oscillations result from neutrino mass. 

The sensitivity of the MINOS near detector neutrinos to a sidereal os- 
cillation signal is shown in Fig. 1 based on a figure from Kostelecky and 
Mewes 8 . 
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Figure 1. The region in this figure labeled "MINOS ND" shows the MINOS sensitivity 
to sidereal oscillation signals resulting from SME. 

2. Data Sample 

The MINOS experiment has been described elsewhere 4 . 

The simulated data set of MINOS near detector neutrinos we use here 
was generated by the standard MINOS Monte Carlo 6 . In total 4.0 x 10 6 
beam spills were generated. In each simulated beam spill we identified the 
charged current (CC) events. These are events with at least one recon- 
structed muon track, as found by the standard MINOS reconstruction. In 
addition, these events must pass the following cuts: (a) Particle ID cut 6 , 
(b) a cut requiring the reconstructed vertex position be > 50 cm from the 
edge of a partial plane or its outline on a full plane, and (c) a cut that 
requires the event be found in the range 1.73 < z < 4.74 m within the 
detector. Both and CC events were accepted. 

Using the total data set, we constructed a histogram of the number of 
CC events/spill. Using this histogram, we simulated 1,000 experiments. We 
generated each experiment by simulating the total number of spills collected 
in the MINOS near detector, 2.78 x 10 6 . For each spill in the experiment 
we picked the number of events found in it from the #CC events/spill 
distribution and we assigned to each a random sidereal time. We then 
binned these events into a sidereal phase histogram spanning 0-1 in sidereal 
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phase, (j> = wqTq. In addition, we binned the protons-on-target (POT) in 
the spill into a second sidereal phase histogram. We then divided the two 
histograms, with the result being a histogram of =#=CC events/POT as a 
function of sidereal phase. Since the #CC events/spill is proportional to 
POT, this final histogram gives the normalized quantity in which we search 
for sidereal variations. Fig. 2 shows a typical Monte Carlo experiment. 
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Figure 2. Left panel: a typical MC experiment; the mean rate is superposed on sidereal 
phase distribution. Right panel: distribution of the fluctuations about the mean rate for 
this experiment; superposed is a Gaussian fit to the fluctuation distribution. 



There are 32 bins in the sidereal phase histograms. Since we are search- 
ing for sidereal variations with an FFT 9 and eq.(l) puts power into fre- 
quencies associated with Fourier terms u>q x n, where n — 1-4, the optimal 
number of bins is given by 2 N with TV = 5. For N — 4, there are too few 
Fourier terms kept in the analysis. 

For each of the 32 sidereal phase bins in a simulated experiment, we 
computed the fractional deviations of the rate in phase bin i, Ri, from the 
mean rate, R, 

Ai?, = ^. (2) 

The distribution of these fluctuations ARi for a typical experiment is shown 
in Fig. 2 with a Gaussian fit superposed. This figure shows that the distri- 
bution is consistent with statistical fluctuations and suggests that residual 
sidereal systematics in the simulated experiments have been minimized. 
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3. Results 

We performed an FFT 9 analysis on each of the 1,000 simulated experiments 
without a sidereal signal. The distribution of the powers in the even and 
odd terms out to frequencies w ffl x n, where n — 1-4, are shown on the left 
and right panels of Fig. 3, respectively. Superposed on these distributions 




FFT Power (cos) FFT Power (sin) 

Figure 3. The results of an FFT analysis of the 1,000 simulated experiments without 
a sidereal signal. The left panel shows the distribution of powers in the even terms and 
the right panel shows the distribution of the odd terms. Only terms out to ojq x n, 
where n = 1-4, are represented in the histogram. Superposed on these distributions are 
Gaussian fits of widths a = 6.3 X 10 -2 for the even terms and cr = 6.6 X 10 — 2 for the 
odd terms. 

are Gaussian fits. We use these distributions as a measure of the signal 
power required for a detection. In this analysis, a detection was defined as 
power in a Fourier expansion term that falls 5<r from the mean as defined 
in Fig. 3. 

We can now set limits on the individual Lorentz and CPT violating 
terms (az,)^ and {cl)^ v making up the coefficients A c , A s , B c , and B s in 
eq.(l) as follows. We first set all (ar J )' 1 and (cr J )'" y terms equal to zero 
except one. Then we increase the magnitude of the nonzero term until 
there is an oscillatory signal with power in a Fourier component that falls 
5(7 from the mean. Fig. 4 shows the 5<r signal that results from the nonzero 
Lorentz violating term 3 (a,L) x = 1.7 x 10~ 19 in both the sidereal phase 
diagram and in the power spectrum. 

4. Future Work 

In the near future we expect to apply these analysis methods to the MINOS 
near detector CC neutrino event data sample. We then plan to extend 
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Figure 4. The 5(7 oscillatory signal resulting from the nonzero Lorentz violating term 
(a^) x . The left panel shows the effect of this signal on the sidereal phase diagram. The 
right panel shows the 5<r deviation in the power in the sin(2LU^T^) term. 



this analysis to include the neutral current (NC) neutrino events. These 
additional events increase the statistics in the data sample by approximately 
50%, thereby improving the significance of the analysis. 
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